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Abstract 

A Dalitz plot analysis of the D° hadronic decays D° -> K°K--k + , D° -» K°K + tt~ and L>° -» 
K K + K~ is presented. This analysis is based on a data sample of 22 fb _1 collected with the BABAR 
detector at the PEP-II asymmetric-energy B Factory at SLAC running on or near the T(4S) 
resonance. The events are selected from continuum e + e _ annihilations using the decay D* + — > 
D°ir + . Preliminary measurements of the branching fractions of the above hadronic decays are 
obtained. Preliminary estimates of fractions and phases for resonant and nonresonant contributions 
to the Dalitz plot are also presented. 
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1 Introduction 



The Dalitz plot analysis is the most complete method of studying the dynamics of three-body 
charm decays. These decays are expected to proceed through intermediate resonant two-body 
modes [||] and experimentally this is the observed pattern with some important exceptions. In the 
case of the decay D + — » K~tt + it + Q, for example, the data can be described with a large 90%) 
nonresonant contribution. Dalitz plot analyses can provide new information on the resonances that 
contribute to observed three-body final states. 

In addition, since the intermediate two-body modes are dominated by light mesons, new in- 
formation on light meson spectroscopy can be obtained. In particular, old puzzles related to the 
parameters and the internal structure of several light mesons can receive new inputs. 

This paper focuses on the study of three-body D° meson decays^] involving a K° (where K° — > 
7r + ir~) such as 

(1) 
(2) 

(3) 

(4) 

K°Ktt. All decays are collectively 



2 The BABAR Detector and Dataset 

The data sample used in this analysis consists of 22 fb _1 recorded with the BABAR detector at 
the SLAC PEP-II storage ring between October 1999 and December 2000. The PEP-II facility 
operates nominally at the T(4S) resonance, providing collisions of 9.0 GeV electrons on 3.1 GeV 
positrons. The data set includes 19.6 fb _1 collected in this configuration (on-resonance) and 2.4 
fb _1 collected below the BB threshold (off-resonance). 

A more complete overview of the BABAR detector can be found elsewhere §. The following is a 
brief description of the components important for this analysis. The interaction point is surrounded 
by a 5-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) filled 
with a gas mixture of helium and isobutane in a 1.5 T superconducting solenoidal magnet. In 
addition to providing precise spatial hits for tracking, the SVT and DCH also measure dE/dx, 
which provides particle identification for low-momentum charged particles. At higher momenta 
(p > 0.7 GeV/c) pions and kaons are identified by Cherenkov radiation observed in the DIRC, a 
detector designed to measure internally reflected Cherenkov light. The typical separation between 
pions and kaons varies from 8a at 2 GeV/c to 2.5a at 4 GeV/c, where a is the experimental 
resolution for the measurement of the Cherenkov angle. 



4 A11 references in this paper to a specific charged state, unless otherwise specified, imply also the charge conjugate 
state. 



£>° -» ITV+Tr-, 
D o ^ K°K-7r + , 

D° — > K°K + ir~, 

D° ->K°K + K-. 

When decays (||) and (J3j) are combined, they are labelled as D' 
referred to as K°h + h~. 
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3 Event Selection and D° Reconstruction 



The decay D* + — > D°ir + is used to distinguish between D° and -D and to reduce background. For 
example, the Cabibbo-favored decays are 

D* + -» J D°vr+ 

-» A' vr + 7r-. 

The charge of the slow vr 1 * 1 from D* decay (referred to as the slow pion) identifies the flavor of the 
D° and K° (apart from a small contribution from doubly-Cabibbo-suppressed decays). The decay 
channels K°K~ir + and K°K + ir~ are produced by different decay diagrams, as shown in Fig. |l|; 
therefore their rates are not expected to be the same. 




Figure 1: Diagrams that contribute to (a) D° — ► K°K 7r + and (b) D° — > K K + ir . 



D° -> K°h+h- candidates are reconstructed from A^ — > 7r + 7r candidates plus two charged 
tracks, each with at least 12 hits in the DCH, and a track with momentum smaller than 0.6 GeV/c 
with at least 6 hits in the SVT. In addition, tracks are required to have transverse momentum 
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Pt > 100 MeV/c and, except for the K° decay pions, to point back to the nominal interaction 
point within 1.5 cm transverse to the beam and 3 cm along the beam direction. 

K° particles are reconstructed by kinematically fitting all pairs of positive and negative tracks. 
Reconstructed K° candidates are then fit to a common vertex with all remaining combinations of 
pairs of positive and negative tracks. Fake K° candidates are removed requiring a flight distance 
of 0.4 cm with respect to the candidate D° vertex. The D candidate is then combined with all 
the slow pion candidates, which are refitted to a common vertex constrained to be located in the 
interaction region. In all fit probability cut at 0.1 % is applied. 

In order to reduce the combinatorial background, D° particles originating from B decays are 
rejected by requiring that the center of mass momentum of the D° candidate is greater than 2.2 
GeV/c. 

Two different particle identification algorithms (PID) are used for kaon identification, both of 
which take advantage of the measurement of dE jdx in the tracking detectors and a measurement of 
the Cherenkov angle in the DIRC. The first employs less stringent requirements and has an efficiency 
above 95% for kaons with momentum less than 3 GeV/c and ir misidentification probability of about 
20%. The second algorithm is more strict and has kaon identification efficiencies of 70% to 90% 
in the same momentum range but with ir misidentification probabilities below 7%. The more 
permissive K identification algorithm is used for the measurement of branching ratios, where more 
uniform acceptance is desirable, and the tighter K algorithm is used for the Dalitz plot analyses, 
where background reduction is more important. 

Each D° sample is characterized by the distributions of two variables, the invariant mass of the 
candidate D°, and the difference between the invariant masses of the D* and D° candidates, 

Am = m(K°h + h--K+) - m(K h + h-), 

where 7i\t is the slow pion. The distribution of Am for those candidates that fall within 2.5 a of the 
D° mass peak is shown in Fig. 0. A strong D* signal is apparent. Fits to these distributions produce 
consistent means and widths for the four channels (a = (323 =L 10) keV/c 2 in the case of decay 
channel K tt + 7t~). The K°h + h~ mass distribution for candidates that fall within ±969 keV/c 2 
(corresponding to 3 a) from the central value of the Am distribution is shown in Fig. |3|. 

4 Efficiency 

The efficiency for the decays in the four D° samples is determined from a sample of Monte Carlo 
events in which each decay mode is generated according to uniform phase space (such that the 
Dalitz plot is uniformly populated). These events are passed through a full detector simulation 
based on GEANT3 § and subjected to the same reconstruction and event selection as the data. The 
distribution of these events in the Dalitz plot after selection is used to determine the reconstruction 
efficiency. Typical Monte Carlo samples used to compute these efficiencies consist of 400 xlO 3 
events. 

This Dalitz distribution is divided into small cells and fit to a third-order polynomial in two 
dimensions. Cells with less than 50 events are ignored in the fit. The resulting x 2 per number 
of degrees of freedom (x 2 /NDF) is typically 1.1. The fitted efficiencies are shown in Fig. |j. The 
average over the Dalitz plot ranges from 7.3% to 8.7%, depending on the decay mode. 
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Figure 2: Am distributions for D° — ► K°h + h~ candidates, for events in which the K°h + h~ 
invariant mass is within 2.5 a of the D° mass. The lines represent the range of Am used to select 
the D° candidates. 



5 Branching Fractions 

Since all four decay channels K h + h~ have similiar topologies, the ratios of branching fractions, 
calculated relative to the K tt + tt^ decay mode, are expected to have relatively small systematic 
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Figure 3: K°h + h mass distributions for events in which Am is within 969 keV/c 2 of the mean 
Am value for signal events. 



uncertainties. These ratios are evaluated as 

where N{(x,y) represents the number of events measured for channel % and ei(x,y) is the corre- 
sponding efficiency in a given Dalitz plot cell (x, y) . 

In order to obtain the yields and measure the relative branching fractions, the K°h + h~ mass 



12 



I/O 

K 7T7T 



K°Ktt 




Figure 4: Dalitz plot efficiency (x 10) for D° -» *r°7r + 7r - , D° -> K^vr and L>° -» K°K + K~ 



distributions are fit with a Gaussian function and a linear background. The number of signal events 
is calculated as the difference between the total number of events within 3a of the D° mass from 
the fit and the integrated linear background function in the same mass range. The results obtained 
with this procedure are summarized in Table [j]. 

Systematic errors take into account the different algorithms used for particle identification and 
the way in which the Dalitz plot efficiency is used to correct the data. For example, the full data 
set is divided into two samples corresponding to different DCH high voltage operating points and 
the branching ratios are computed separately for the two samples. Differences in the branching 
ratios are used as an estimate of systematic uncertainties. The contributions to the systematic 
errors are summarized in Table The total error is obtained by adding in quadrature the single 
contributions. 

The method for measuring the branching ratios is checked with a different fully inclusive e + e~ — > 
cc Monte Carlo sample in which the D° mesons decay according to the PDG || branching fractions. 
The Monte Carlo events are subjected to the same reconstruction, event selection and analysis as 
the real data. The size of this Monte Carlo sample is comparable to that of the data sample. The 
results are found to be statistically consistent with the branching ratios assumed in the Monte 
Carlo generation. 
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Table 1: Results from the fits to the K°h + h~ mass distributions. The fits are performed using a 
linear background and a Gaussian function for the signal. 



Channel 




mass (MeV/c 2 ) 


a (MeV/c 2 ) 


signal events 


D° - 




IT 


1863.9 ± 0.5 


6.29 ± 0.05 


15 279 ± 129 


D° - 


K°K' 


-7T+ 


1863.9 ± 0.2 


5.1 ± 0.2 


1 165 ± 40 


D°- 


,K°K- 




1864.3 ± 0.2 


5.0 ± 0.2 


805 ± 35 


D°- 






1864.7 ± 0.8 


3.6 ± 0.1 


2 109 ± 48 



Table 2: Contributions to the systematic error on the ratio of branching fractions relative to 

7T + 7f — (%) from different sources. 

Channel PID Efficiency correction DCH voltage Total 





-> K 1) K-tt+ 


0.31 


0.14 


0.45 


0.56 


D° - 


-»• K°K+Tr- 


0.18 


0.04 


0.37 


0.41 


D° - 


-> K°K+K- 


0.19 


0.09 


0.17 


0.27 



The measured preliminary ratios of branching fractions are shown in Table |3| and are compared 
with those measured by other experiments. 



Table 3: Ratio of branching fractions relative to D c 
to K°K~tt + from if V +tT 



TtO 4- 



(%). Note that ARGUS did not 



se parate A 
Channel 



CLEO § ARGUS § 1X39T 



BABAR 



(T 



D 



D° -> K°K+K- 



8.32 ± 0.29 (stat) ± 0.56(syst) 
5.68 ± 0.25(stat) ± 0.41(syst) 
16.30 ± 0.37(stat) ± 0.27(syst) 



10.8 ± 1.9 16.0 ± 3.0 12.4 ± 4.7 
9.8 ± 2.0 7.8 ± 3.9 

17.0 ± 2.2 



6 Dalitz Plot Analysis Method 

An unbinned maximum likelihood fit is performed for the decay modes D° — ► K°K~tt + , D° — » 
K K + TT~ and D° -> K K+K~, in order to use the distribution of events in the Dalitz plot to 
determine the relative amplitudes and phases of intermediate resonant and nonresonant states. 

Following the method used by ARGUS || and CLEO [ 10 1 , the likelihood function is written in 
the following way: 



L=p c{m) >:,,'■>'■} a > a j , a-/?) (5) 

/ "52i,j CiC*AiA*e(ml,rn^)dmldm^ J e(ml,m^)dmldm^ 

In this expression, (3 represents the fraction of signal obtained from the fit to the mass spectrum and 
e(m^,mH) is the Dalitz plot efficiency. G(m) is a Gaussian function describing the D° lineshape 
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normalized within the 2a cut used to perform the Dalitz plot analysis. It is assumed that the 
background events, described by the second term in Eq. ||, uniformly populate the Dalitz plot. 
This assumption is verified by examining events in the D° side bands. The output from the fit is 
the set of complex coefficients Cj. 

In Eq. 5, the integrals are computed using Monte Carlo events taking into account the Dalitz 
plot efficiency. The branching fraction for the resonant or nonresonant contribution i is defined by 
the following expression: 

|q| 2 J \ Ai\ 2 dmldm,y 



f. = 



2 ' 



The fractions fi do not necessarily add up to 1 because of interference between amplitudes. The 
errors in the fractions are evaluated by propagating the full covariance matrix obtained by the fit. 

The phase of each amplitude is measured with respect to the mode with the largest amplitude. 
The amplitudes Ai are each represented by the product of a complex Breit-Wigner function and 
an angular function (the Zemach tensors [pd| ): 

A = BW(m) x T(f2). 



The Breit-Wigner function includes Blatt-Weisskopf form factors [12]. The /o(980) and ao(980) 
resonances is described using coupled-channel Breit-Wigner functions with parameters taken from 
the CERN/WA76 Jl^] and the Crystal Barrel []l4|] experiments, respectively. The parameters of the 
0(1020) meson are extracted from the data (r = (4.3 ± 0.3) MeV/c 2 ) since the apparent width is 
affected by the experimental resolution. The resonance parameters of the ^(1430) are m = 1.435 
GeV/c 2 and T=279 GeV/c 2 from a recent reanalysis of elastic Kir scattering data from the LASS 
experiment [[15|| . The nonresonant contribution (N.R.) is represented by a constant term with a 
free phase. 

The fit quality is evaluated in the following way. The Dalitz plot is re-binned into N ce u s cells 
grouping together bins with small event yields. Then a Monte Carlo sample is generated with 
parameters from the fit to data. This distribution is used to evaluate a x 2 f° r t ne Dalitz plot. The 
X 2 /N cells is displayed in Table |I] for the fit results shown in Tables ||, ||, and 0, together with the 
event yield and purity for each of the three channels used in the Dalitz plot analyses. 



Table 4: Event yields, purity and X 2 l^ceiis f° r the three channels in the Dalitz plot analyses. 



D° Decay mode 


Events 


Purity (%) 


X 2 /N cells 


D u K°K-ir+ 
D° -> K°K+7r- 
D° -> K°K + K~ 


1008 

659 

1957 


95.5 ± 0.4 
95.5 ± 0.4 
97.5 ± 0.2 


46/44 
25/29 
98/59 



Systematic errors on the fitted fractions are evaluated by making different assumptions in the 
fits — for example, we have performed fits with the Blatt-Weisskopf terms set to 1, and with 
uniform efficiency across the Dalitz plot. 
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7 D° —*■ K°K-n + Dalitz Plot Analysis 

The Dalitz plot for D° — > i^°X~7r + candidates is shown in Fig. [5| and its projections are shown in 
Fig. ^. The presence of a strong K* + (892) resonance can be observed. Table [5] shows the list of 
intermediate final states and the fitted fractions for this decay channel. 



2 



> 1 5 

CD 



CM 



0.5 



^ V ; • ... 



J L 



0.5 1 1.5 2 

m 2 (K° n + ) GeV 2 /c" 



Figure 5: Dalitz plot of D° -» K^K~ix 



The results from the fit, shown in Table [El, confirm that the channel is dominated by K* 
resonances in the charged mode. Contributions from neutral K* as well as contributions from ao 
resonances are small or consistent with zero. The nonresonant contribution appears to be negligible. 

Recently, in order to have a better description of the Dalitz plot of D + — > iT - 7r + 7r + , experiment 
E791 |Q] introduced in the analysis the old postulated Kir S-wave resonance k fl6|| . The best 
parameters found were M K = (797 ± 19 ± 42) MeV/c 2 and T K = (410 ± 43 ± 85) MeV/c 2 . Including 
a contribution from n + K~ in the present analysis results in an increase of 21ogL of 7.8 for 2 more 
parameters. The resulting fraction for this contribution is (15 ± 12)%. The large error does not 
allow us to confirm the existence of this final state. 
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Table 5: Preliminary results from the Dalitz plot analysis of D° — > K°K it 



Final state 


Fraction 


(%) 


Phase (degrees) 


K* °(U30)K° 


4.8 ± 1.4 


± 1.6 


52 ± 27 


Kl°(892)K° 


0.8 ± 0.5 


± 0.1 


175 ± 22 


Kf (1680)K° 


6.9 ± 1.2 


± 1.0 


-169 ± 16 


Kf(U30)K° 


2.0 ± 0.6 


± 0.1 


51 ± 18 


K* + {U30)K- 


13.3 ± 3.5 


± 3.9 


-41 ± 25 


K* + (892)K- 


63.6 ± 5.1 


± 2.6 





Kl + (1680)K~ 


15.6 ± 3.0 


± 1.4 


-178 ± 10 




13.8 ± 2.6 


± 7.9 


-52 ± 7 


a (980)vr+ 


2.9 ± 2.3 


± 0.7 


-100 ± 13 


a (1450)vr + 


3.1 ± 1.9 


± 0.9 


31 ± 16 


a 2 (1310)vr+ 


0.7 ± 0.4 


± 0.1 


-149 ± 27 


N.R. 


2.3 ± 0.5 


± 5.6 


-136 ± 23 


Sum 


130 ± 


8 






12 3' 1 2 1 



m 2 (K° K") m 2 (K° tt + ) m 2 (K~ tt + ) 

Figure 6: Dalitz plot projections for D° — > K K~ir + . The data are represented with error bars; 
the histogram is a projection of the fit described in the text. 

8 D° — ► K°K + n- Dalitz Plot Analysis 

The Dalitz plot for D° — > K°K + tt~ candidates is shown in Fig. [?] and its projections are shown 
in Fig. ||. The presence of the K* + (892) resonance is evident together with a small K*°(892) 
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Figure 7: Dalitz plot of L>° -> K°K + ir 



contribution with asymmetric lobes suggesting interference with other final states. A shoulder near 
the threshold of the K K+ mass distribution suggests the presence of the tig" (980) resonance. This 
D° decay channel is fit using the final states shown in Table ^, which also displays the fit fractions 
and phases. The resulting x 2 f° r this fit is given in Table ||. 

This fit reveals a rather complex structure with large uncertainties due to the limited sample 
size. A sizeable nonresonant contribution is predicted by the fit, in contrast to the D° — » K°K~tt + 
channel discussed in the last section. 

Dropping the nonresonant contribution results in a decrease of 21ogL of 9 units for 2 less 
parameters with the same x 2 value. While all other fractions have little changes, in this fit the 
contribution from K*~(U30)K + increases from (2 ± 8)% to (26 ± 12)% and that of £^(980)7^ 
decreases from (15 ± 13)% to (5 ± 6)%. The sum of the fractions for this fit is (114 ± 23)%. 



9 D° —*■ K°K + K~ Dalitz Plot Analysis 

The Dalitz plot for D° —* K°K + K~ candidates is shown in Fig. M and its projections are shown 



in Fig. 10. A strong ^(1020) signal interfering with a threshold scalar meson can be clearly seen. 
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Figure 8: Dalitz plot projections for D° — > K K + ir . The data are represented with error bars; 
the histogram is a projection of the fit described in the text. 



Both /o(980) and a[j(980) resonances can be present near threshold. An accumulation of events 
due to a charged (980) can be observed on the lower edge of the Dalitz plot. Table |7| shows the 
list of the final states and the fitted fractions and phases for this D° decay channel. 

Most of the uncertainty in this channel is due to the poorly defined parameters of ao(980) 
and /o(980). Their properties are taken from measurements performed by other experiments. 
Since these states lie below KK threshold, their properties cannot be measured in this D° decay 
mode only. Their line shape would be best determined from a coupled-channel analysis of K^K^, 
KgK®, K + K~, rjir and tttt decays. A scan is made of the likelihood function with respect to the 
ratio g K x/gr)n for ao(980) leaving the other parameters fixed. The likelihood has a maximum for 
g K xl g-q-K ~l-3, to be compared with the Crystal Barrel result of 1.01 ± 0.07. 

Similarly, a scan is performed as a function of gx/g-ir for /o(980). The likelihood has a maximum 
for gx/g-K ~1.25, to be compared with the WA76 result of 2.0 ± 0.6. The variation of the fractions 
due to the change of these parameters is taken into account in the evaluation of the systematic 
errors. 

The doubly-Cabibbo-suppressed contribution D° — > K°K + K~ has also been included in the 
fit. The presence of such a contribution should show up as an Oq (980) signal in the K°K~ effective 
mass. The insertion of this new term results in an increase of 2 log L by 10 units for 2 additional 
parameters and the resulting \ 2 is included in Table ||. However, with the present statistics, its 
contribution of (3.2 ± 1.9)% is consistent with zero. 

The poor fit quality for this channel, indicated by the X 2 /N ce u s of 98/59 in Table ||, comes 
mostly from the (ft mass region. Effects of (ft mass resolution comparable to its natural width 
combined with uncertainties in the ao and /o line shapes are not included in the current fit model, 
and do not correctly model the interference in this region. Future analyses with larger data samples 
should attempt to address this problem with a more detailed model. 
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Table 6: Preliminary results from the Dalitz plot analysis of D° — > K K + ir . 
Final state Fraction (%) Phase (degrees) 



^°(1430)K° 


26.0 


± 


16.1 ± 3.3 


-38 ± 22 


iff (892)1?° 


2.8 


± 


1.4 ± 0.5 


-126 ± 19 




15.2 


± 


11.9 ± 0.5 


161 ± 9 


Kf(U30)K° 


1.7 


± 


2.5 ± 0.2 


53 ± 38 


K*-{U30)K+ 


2.4 


± 


8.2 ± 1.0 


-142 ± 115 


K*-(892)K+ 


35.6 


± 


7.7 ± 2.3 





(1G80)K+ 


5.1 


± 


5.7 ± 1.1 


124 ± 27 




1.0 


± 


1.0 ± 0.2 


-26 ± 38 


a+ (980)tt- 


15.1 


± 


12.5 ± 0.6 


-160 ± 42 


a+(1450)vr- 


2.2 


± 


2.7 ± 1.2 


148 ± 25 


N.R. 


36.6 


± 


25.8 ± 2.7 


-172 ± 13 



Sum 144 ± 37 



Table 7: Preliminary results from the Dalitz plot analysis of D° — > K°K + K 

Final state Fraction (%) Phase (degrees) 

K°(f) 45.4 ± 1.6 ± 1.0 

K°a° (980) 60.9 ± 7.5 ± 13.3 109 ± 5 

K°/o(980) 12.2 ± 3.1 ± 8.6 -161 ± 14 

a (980)+K- 34.3 ± 3.2 ± 6.8 -53 ± 4 

a (980)-K+ 3.2 ± 1.9 ± 0.5 -13 ± 15 

N.R. 0.4 ± 0.3 ± 0.8 40 ± 44 

Sum 156 ± 9 
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Figure 9: Dalitz plot of L>° -> K K+R-. 



10 Summary 



Dalitz plot analyses are performed for the hadronic decays D° — > K°K 7r + , D° — » i^°X + 7r and 
Z)° — > i<C K + K~ . The fractions and relative phases for intermediate resonant states are extracted. 
The following preliminary ratios of branching fractions are measured: 

(8.32 ± 0.29(stat) ± 0.56(syst)) x 1(T 2 , 

(5.68 ± Q.25(stat) ± OAl(syst)) x 1(T 2 , 
(16.30 ± 0.37 (stat) ± 0.27(syst)) x 1(T 2 . 
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Figure 10: Dalitz plot projections for D° — > K°K + K . The data are represented with error bars; 
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